Astronomy & Astrophysics manuscript no. kurk gmassV zspike 


© ESO 2009 


June 24, 2009 





GMASS ultradeep spectroscopy of galaxies at z - 2"^ 

V. Witnessing the assembly at z = 1.6 of a galaxy cluster 

J. Kurk^'^, A. Cimatti^ G. Zamorani"^, C. Halliday^, M. Mignoli"^, L. Pozzetti^, E. Daddi^ P. Rosati^, M. Dickinson^, 
M. Bolzonella"^, P. Cassata^, A. Renzini^, A. Eranceschini^^, G. Rodighiero^^, and S. Berta^^ 

^ Max-Planck-Institut fiir Astronomie, Konigstuhl 17, D-69117, Heidelberg 

e-mail: kurk@mpia.de 
^ INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 1-50125, Firenze 



Universita di Bologna, Dipartimento di Astronomia, Via Ranzani 1, 1-40127, Bologna 
^ INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1, 1-40127, Bologna 

^ CEA, Laboratoire AIM - CNRS - Universite Paris Diderot, Irfu/SAp, Orme des Merisiers, F-91 191 Gif-sur-Yvette, France 
^ European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85748, Garching bei Munchen 
^ NOAO-Tucson, 950 North Cherry Avenue, Tucson, AZ 85719, USA 

^ Department of Astronomy, University of Massachusetts, 710 North Pleasant Street, Amherst, MA 01003, USA 

^ INAF-Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 1-35122, Padova 
Universita di Padova, Dipartimento di Astronomia, Vicolo dell'Osservatorio 2, 1-35122, Padova 
Max-Planck-Institut fiir extraterrestrische Physik, GiessenbachstraBe, D-85748, Garching bei Munchen 
Received ; accepted 

ABSTRACT 

Context. Clusters of galaxies represent important laboratories for studying galaxy evolution and formation. Well established and 
relaxed clusters are known below z < 1.4, as well as, clusters in formation found close to radio galaxies at z > 2, but the in-between 
redshift range, during which clusters are expected to undergo significant changes, is almost unexplored. 

Aims. By studying a galaxy overdensity in redshift and angular distribution at z = 1.6, uncovered in the Galaxy Mass Assembly 
ultra-deep Spectroscopic Survey (GMASS), we provide insight into the evolution of cluster galaxies at high redshift. 
Methods. We present a study of the significance of the galaxy overdensity at z = 1.6, CI 0332-2742, its velocity dispersion, and X-ray 
emission. We reveal colour bimodality for cluster members and compare the properties of members of CI 0332-2742 with galaxies 
outside the overdensity. 

Results. From the redshifts of the 42 overdensity members, we find a velocity dispersion of 500 kms~^ We conservatively estimate 
the overdensity in redshift space for the spike at z = 1.6 in the GMASS to be 8.3+1.5. A map of the surface density of galaxies 
at z = 1.6 in the GMASS field shows that its structure is irregular with several filaments and local overdensities. The diff'erences 
in the physical properties of CI 0332-2742 member and field galaxies agree with the latest hierarchical galaxy formation models: 
for overdensity members, the star formation rate (SFR), and specific SFR, is approximately 50% lower than for the field galaxies; 
overdensity galaxies are twice the age, on average, of field galaxies; and there is a higher proportion of both massive (M > 10^^ ^M©), 
and early-type galaxies, inside CI 0332-2742 than in the field. Among the 42 members, seven have spectra consistent with being 
passively evolving, massive galaxies. These are all located within an area where the number density of z = 1.6 galaxies is highest. 
In a z - / colour-magnitude diagram, the photometric data of these early-type galaxies are in close agreement with a theoretical red 
sequence of a galaxy cluster at redshift z = 1.6, which formed most of its stars in a short burst of star formation at z ~ 3. 
Conclusions. We conclude that the redshift spike at z = 1.6 in the GMASS field represents a sheet-like structure in the cosmic web 
and the area with the highest surface density within this structure, containing already seven passively evolving galaxies, will evolve 
into a cluster of galaxies at a later time. 

Key words. Galaxies: distances and redshifts - Galaxies: evolution - Galaxies: formation - Galaxies: clusters: general - Galaxies: 
high-redshift 



1. Introduction 
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Clusters of galaxies are the most massive structures in the uni- 
verse. Because the evolution of structure in the early universe 
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depends on mass, galaxy clusters play an important role in cos- 
mology and can help to constrain cosmological parameters (e.g., 
iRapetti et al. 2008). Clusters are also powerful tools for studying 
galaxy evolution. 

It is well known that environment plays an important role 
in galaxy evolution. Correlations between the environment of 
galaxies, e.g., in terms of galaxy density, and several galaxy 
properties, such as the age and metallicity of their stellar popula- 
tions, star formation rates, and morphologies have been studied, 
in individual galaxy clusters and also as a function of redshift. 
Early observations of cluster galaxies o ut to z = 0.5 by, for ex- 
ample, iButcher & Oemled (Il978l 11984 showed little evidence 
of evolution in the colours of the red cluster population, but 
detected an increase in the fraction of blue galaxies, which ap- 
peared to avoid the cluster centre. Although this Butcher- Oemler 
effect has been disputed, the general conclusion still st ands (e.g., 
Ellingso n etaL''2001: 'Ta naka etal .' '2005: "Andreon et al. 200^ 
Tran et al. 2007 ; Fassbender et al.i „2008). Balogh et al.i (.20071) 
find that, for a given stellar mass and redshift, fewer galaxies 
in groups show strong [Oil] emission than in field galaxies. 
The evolution of star formation rate appears to dep end on both 
galaxy mass an d galaxy envi ronment (Bal dry et al.l 12006) . The 
observations by iBaldry et al.l of a large sample of galaxies at 
z < 0.1 were compared directly with the se mi-a nalytic mod- 
els of galaxy formation by lBower et al.l (1200 6) and'C roton et al.l 
(E006), showing closer agreement with the former. This colour- 
density r elation has also been confirmed at higher redshifts, up to 
z ~ 1.5 (ICucciati et al.ll2006l: ICooper et al.ll2QQ7l). At these high 
redshifts, few massive systems are known, as iPoggianti et al.l 
(2008) note. In addition, distant cluster surveys have studied pri- 
marily massive clusters. Poggianti et al., therefore, employ data 
from the ESO Distant Cluster Survey to investigate the relation 
between star formation activity, morphology and local galaxy 
density in galaxy groups and clusters at z = 0.4-0.8. They de- 
tect proportionally fewer galaxies with ongoing star formation 
in regions of higher projected density in both distant and nearby 
clusters, but no evidence that the star formation properties of a 
given Hubble type or galaxy mass depend on local environment. 

A striking feature in the spectrum of an old stellar popu- 
lation is the 4000 A break, above which the light of old (red) 
stars dominate over the rare young stars. This feature causes the 
occurrence of a sequence of red (passive) galaxies in a (suit- 
able) colour-magnitude diagram of galaxy clusters. The tight 
sequence observed in both nearby and distant clusters (up to 
z = 1.24) suggests that the stars in these gal axies formed at very 
high redshift in a relatively short time (e.g., Stanford et al.ll 19981 
iBlakeslee et al.l l2003). An analysis of stellar population proper- 
ties and morphologies of red-sequence galaxies in clusters and 
groups from z ~ 0.75 to z ~ 0.45 by Sanch ez-Blazquez et al.l 
(I2009h shows that the rate at which red-sequence galaxies evolve 
depends on their mass, as less massive galaxies show evidence 
for a more extended star formation history than more massive 
galaxies do. One therefore expects that, at higher redshift, the 
scatter in colour increases as the epoch when the stellar pop- 
ulations were formed is approached. Indeed, changes in the 
red sequence in galaxy structures have been observed between 
low redsh ift and redshifts z = 2.2-3.1 (Kodamaet al. 20071: 
IZirmetal.ll2008h . in the sense that the faint end of the red se- 
quence becomes less populated, although a depe ndence on clus- 
ter richness may also play a role (lAndreonll2008h . 

Distant groups and clusters of galaxies, therefore, provide 
important environments for the study of galaxy evolution, es- 
pecially if both passive and actively star-forming galaxy popu- 
lations can be studied in detail. For the overdensities of galax- 



ies known at high redshift, this is often problematic because 
their cluster members are selected on the basis of either their 
star formation activity (blue members, e.g., Kurketal. 2000; 
IPentericci et al .120001: 1 Venemanset al . 2002; Steidel et al. 2005|), 
or lack of it (red members, e. g., iBest et al.„ 2003: >Kodama et al.l 
120071: [McCarthy et"an l2007l) . The most well studied high- 
redshift systems confirmed to be clusters on the basis of their 
X-ray emission, are found between z = l.O and z = 1.4, red- 
shifts for which spectroscopic data of both the red and blue 
populations can be practically acqu ired (e.g., at z = 1.26, 
iDemarco etaP 120071: Mei et al.ll2006h . These clusters were dis- 
covered either by their extended X-ray emission due to hot clus- 
ter gas (RX J0848.9-h4452 at z = 1.26, XMMU J 2235.3-2557 
at z = 1.39, XMMXCS J 2215.9-1738 at z = 1.45, iRosati et al.1 
119991: iMullis et al.ll2005l: [Stanford et al.ll2006[ respectively), or 
by studying near-infrared (NIR) photom etry of red cluster me m- 
bers (ISCS J143809-F341419 at z = 1 .41 [Stanford etaDl2005h . 

We present a galaxy overdensity at z = 1.61, possibly a clus- 
ter under assembly, containing both red and blue galaxy pop- 
ulations. In Sect. El we introduce the overdensity as described 
in other papers and describe the data set that we used to gather 
additional data. We also discuss the selection bias that may be 
present in our data. In Sect. [3j we describe the properties of the 
overdensity, including its velocity dispersion, an estimate of its 
mass, and the presence of a red sequence of galaxies. In Sect.lH 
the properties of the individual galaxies in the overdensity are 
described, in relation also to the local density, and in comparison 
with galaxies in a field sample. The overdensity and its galaxies 
are discussed in the context of the literature in Sect. [51 In Sect.[6l 
the paper is summarized. 

The following cosmological parameters are used throughout 
this paper: Hq = 70 km s"^ Mpc-\ = 0.3, Qa = 0.7, Qk = 0. 
Magnitudes are quoted in the AB system. 

2. The data 

2.1. An overdensity at z =1.6 

The overdensity at z = 1.6 was detected by means of its clus- 
tering in the redshift space of the deep spectroscopic GMASS 
survey (Kurk et al., in preparation). We refer to this overden- 
sity as CI 0332-2742, which corresponds to the coordinates of 
the galaxy at the centre of the number density isosurfaces (see 
below). The GMASS survey was carried out in the field of 
GOODS-S, w here several overd ensities at various redshifts have 
been detected. [GiUi et al.[ ([2003[) described five spikes in the X- 
ray source redshift distribution of this field; the second high- 
est in terms of redshift occurs at z = 1.62 and contains five 
sources. The Poissonian probability of observing such a high 
number of sources, given the background count rate, is only 
4 X 10"^, whic h implies that t his aggregation represents a sig- 
nificant group. [Vanzella et al.[ ([2006) also described spikes in 
the redshift distribution using data from the ESO/GOODS spec- 
troscopy program of faint galaxies in the GOODS-S field. They 
reported that 20 galaxies have redshifts of z ~ 1.61 and eight 
other galaxies had redshift data from other surveys; all of these 
galaxies are distributed spatially in an apparently non-uniform 
way. The probability of detecting such a significant peak in 
redshift space by chance is less than 7 x 10"^, according to 
Monte Carlo simulations of the smoothed redshift distribution 
of the 501 confirmed redshifts, after the second stage of the 
ESO/GOODS spectroscopy program. 

Three old, fully-assembled, massive (> 10^ ^M©) spheroidal 
galaxies within the spike at z = 1.61 were detected in the 
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K2Q s urvey (ICimatti et al.ll2002h and described in 
20041). n another article based on GMASS results, 
2008) investigated the properties of thirteen old, passive galax- 
ies at z > 1.4, seven of which are at z = 1.6. These galaxies have 
stellar-mass surface-densities that are ~ 1 dex higher than local 
spheroids of similar mass. 

An overdensity at z = 1.6 was also discerned in photo- 
metric redshift space by ICastell ano et al.l (12007); these authors 
applied their own (2-F 1 )D cluster finding algorithm (see also 
iTrevese et aDl2007l) to the GOODS-South Field and discovered 
a localized overdensity, embedded in a larger scale overdensity 
of galaxies at z = 1.6. Their method was based on an esti- 
mate of three-dimensional densities and considered simultane- 
ously angular positions and photometric redshifts (in this cas e 
from the GOODS-MUSIC catalogue by 'Grazianetal] l2006h . 
From their analysis, several major peaks at different redshifts 
emerged, the peak at z ~ 1.6 being the most significant at z > 1 
in the GOODS-South field. Within this large-scale overdensity 
(for which additional spectroscopic evidence was presented by 
ICimatti etaLlEOQl I2QQ41: iGilfi et al.ll2003l: IVanzella et al]l2005[ 
2006), they isolated a compact, higher density peak, centred ap- 
proximately on RA=03^32^29.98^ DEC=-27°42'35.99'^ which 
was identified to be a cluster of an approximate total extension 
of 3 X 3 Mpc (comoving), centred on a large star-forming, re - 
markably irregular galaxy (see Fig. 2 in ICastellano et al.ll2007l) . 
In this paper, we define a circular high-density region, centred 
on the same coordinates, of radius 1 Mpc (physical). The rea- 
sons for this choice are explained in Sect. 13.11 On the basis of 
the irregular morphology of the galaxy overdensity, its low total 
X-r ay emission, and its estimated mass of(1.4-4.1)xl0^'^ M©, 
ICastellan o et al. concluded that this group/poor cluster had not 
yet reached its virial equilibrium. 



2.2. The GMASS sample 

GMASS {"Galaxy Mass Assembly ultra-deep Spectroscopic 
Survey''^ is a project based on an ESO VLT Large Program 
(173.A-0687, P.I. A. Cimatti). The GMASS project and sample 
are described in detail by J. D. Kurk et al. (2009, in prepara- 
tion, hereafter Paper VI), and we recall only its main features. 
The GMASS sample was selected from a reg ion of 6.8 x 6.8 
arcm in^, located in the GOODS-South fielcn (iGiavalisco et al.l 
12004. Dickinson et al., in preparation); all sources detected in 
the publicly available 4.5 yum image, acquired using IRAC on 
the Spitzer Space Telescope, were extracted to a limiting mag- 
nitude of m4.5 < 23.0 (2.3yuJy), which provided a selection 
related closely to stellar mass. The GMASS sample includes 
1277 unblended objects to m4.5 < 23.0, which have photometry 
from t he NUV to MIR and SED fits obtained using Maraston 
(l2005l M05) templates providing stellar masses, star forma- 
tion rates (SFRs), a nd other galaxy properties (see Sec. |4] and 
iPozzetti et aDl2007l) . 

The GMASS ESO VLT-FF0RS2 optical spectroscopy was 
focused on galaxies pre-selected using a cut in photometric 
redshift of Zphot > 1.4 and two cuts in the optical magnitudes 
(5 < 26.5, /< 26.5). This selection produced 221 spectroscopic 
targets, 174 of which were actually observed (called the spectro- 
scopic sample, hereafter). The integration times were extremely 
long (up to 32 hours per mask and some targets were included in 
multiple masks), and the utility of the spectroscopic data was 
optimized by acquiring spectra at either blue (4000-6000 A) 



^ http : //www. arcetri .astro . it/~ cimatti/gmas s/gmas s . html 
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or red (6000-10000 A) wavelengths, depending on the colours 
and photometric SEDs of the targets, called the blue and red 
masks, respectively, hereafter. Despite the faintness of the tar- 
gets, the GMASS spectroscopy program was very successful: we 
were able to determine redshifts for 150 out of the 174 targeted 
sources (86%), 132 of which turned out to have Zspec ^ 

1.4. This 

implies that fractions of >76% and 88% of all targeted sources 
and those with successfully determined redshifts, respectively, 
are at Zspec > 1-4. In addition to the 18 objects from the spec- 
troscopic sample with Zspec < 1-4, filler objects used to complete 
the masks delivered another 27 redshifts of Zspec < 1 .4 and three 
of Zspec > 1.4. 

Figure [T] shows the distribution of all known spectroscopic 
redshifts (including those of other surveys in GOODS-S, see 
Paper VI) of galaxies within the GMASS catalogue (stars are 
omitted), up to z = 2.9. Several peaks in the distribution are 
prominent, the peak at z = 1.6 being the most significant at 
z > 1.4. Although this peak was previously known to exist (see 
Sect. 12.11 for references), the availability of 136 spectroscopic 
redshifts at z > 1.4 from GMASS allows a far more careful 
study than possible before. The redshift peak at z = 1.6 may in- 
dicate the presence of a galaxy cluster or a more extended galaxy 
overdensity (such as a sheet). Because the dynamical state and 
richness of CI 0332-2742 is a priori unknown, we refer to the 
galaxy aggregation as a galaxy group or cluster in formation and 
to the galaxies considered to be part of the spike in the redshift 
distribution as members or member galaxies. 



2.3. Assessment of the spectroscopic selection bias 

Our criteria for the selection of the spectroscopic sample may 
create a selection bias, which should have a minor influence on 
the redshift distribution in the range 1.4 < z < 1.8 considered 
here. A bias introduced by the success of the redshift determina- 
tion, could be more important. At wavelengths for which [O II] 
is detected, in galaxies at z ~ 1.6, several sky lines of a range of 
intensities, in the observed frame, are present. These lines can, in 
principle, have a serious effect on the measurement of emission- 
line wavelengths, and spurious peaks in the redshift distribution 
can be produced. We find, however, that this eff'ect has a minor 
influence on our redshift determination. 

In Sect. 13.21 we show that the 42 galaxies in the redshift 
range 1.600 < Zspec < 1.622 are likely to be members of 
the galaxy overdensity. Of these, 27 have redshifts determined 
only by GMASS, while ten have redshifts determ ined only by 
ESO/GOODS spectroscopy (VanzeU a et aDl2006h . The remain- 
ing five have redshifts determined b y GMASS spec troscopy, but 
four of these were also observed by I VanzeUa et alJ (two with an 
uncertain redshift, one of them being inconsistent with the more 
accurate GMASS redshift), and one was part of the K20 survey 
(but with an uncertain and inconsistent redshift). We examined 
all 42 spectra in detail and assessed whether redshift determina- 
tion would have been possible without the detection of the [O II] 
line, which was the case for 20 (i.e., 48%) of the galaxy spectra 
with strong UV absorption lines. For eight of these spectra, the 
region about A = 9690 A, which corresponds to the wavelength 
of the [O II] line at z = 1.61, was not observed at all, while for an 
additional three the region was observed, but the [O II] line was 
not detected. A wavelength-dependent selection function is not 
an issue for these galaxy spectra, because the main UV absorp- 
tion lines are observed in an observed-frame wavelength range 
for which no strong sky lines are detected. 
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Fig. 1. Redshift distribution of spectroscopic redshifts of ob- 
jects in the GMASS field displayed with bin sizes of Az = 0.02 
(main panel) and 0.002 (insets). Redshifts determined within the 
GMASS survey are indicated by the grey histogram. Within the 
upper inset, the sky background relevant to [O II] detection is 
indicated by a dashed line. Within the lower inset, which has 
the same bin size but is zoomed in on a smaller redshift range, 
including the spike only, a curve with two Gaussian functions 
fitted to the histogram is overlayed. In addition, a dashed his- 
togram showing the 21 galaxies within the high density region 
(see text) is included. The solid arrow indicates the redshift of 
the brightest member galaxy, while the dashed arrow indicates 
the redshift of the central (see text) galaxy. 



As a second check, we assessed precisely the sky line sit- 
uation about 9690 A. The result is shown in the upper inset of 
Fig. [TJ where we have translated the wavelengths of the sky 
spectrum to the redshift relevant to [O II] detection. We note that 
the region corresponding to 1.600 < z < 1.616 is populated by 
weak sky lines. In contrast, the wavelength range corresponding 
to 1.570 < z < 1.600, for which not a single redshift was deter- 
mined, has far lower sky line contamination. In addition, we note 
that in the redshift range 1.616 < z < 1.622, corresponding to an 
observed wavelength range for which no sky lines are present, 
only one galaxy is detected. 

We conclude that the observed-frame wavelength range af- 
fected by mild sky-line contamination is significantly wider than 
the range encompassing the overdensity redshift peak. Sky lines 
have, therefore, a minor influence, if any, on the measurement of 
redshift for the twelve members with GMASS spectroscopy, for 
which redshift was measured using [Oil], and the ten members 
for which redshifts were measured by the ESO/GOODS survey. 

We finally note that the redshift peak is also present in the 
photometric redshift distribution (see Fig. 1 in paper VI), which 
is subject to diff'erent selection biases and therefore unlikely to 
generate a peak at the same redshift by chance. 



3. Overdensity properties 

3.1. Angular distribution 

Figure [2] shows the geometry of the GMASS field, indicat- 
ing objects with known spectroscopic redshifts by small dots, 
galaxy members (i.e., objects with 1.600 < Zspec < 1.622) by 
crosses or diamonds, and photometric redshifts in the range 
1.43 < Zphot < 1-77 by large dots. There is a distinct con- 
centration of galaxy members towards the northern edge of the 
GMASS field. In addition, among the galaxy members, the red, 
early, and intermediate-type galaxies are most strongly concen- 
trated towards the north. These types were determined from the 
spectra of the galaxies on the following basis: class 1 (early- 
type) have absorption lines only, class 2 (late-type) have strong 
emission lines, and class 1.5 (intermediate-type) have absorp- 
tion lines but some emission lines (most commonly [Oil]) are 
also present. 

The brightest confirmed galaxy member (GMASS 2148 at 
z = 1.609, close to the mean redshift of the spike) is the cen- 
tral galaxy within a striking linear triplet of early-type galax- 
ies, having a length of 20'' (or 169 kpc) in the north-east of 
the GMASS field (see Fig. [3]). The other two members of the 
triplet (GMASS 2196 and 2543) have redshifts z = 1.614 and 
1.612 (respectively) and the largest diff'erence in redshift corre- 
sponds to ~ 575 km s"^. The crossing time of a galaxy in this 
triplet is ~ 3 X 10^ yr, indicating that these systems are possibly 
in the process of merging, in which case they would form one 
central massive galaxy with a mass of 3 x lO^^M©. In a colour- 
magnitude diagram of this field, GMASS 2196 and 2543 are also 
located very close to GMASS 2148, the latter forming the tip of 
the red sequence (see Sect. l3.5b . It is remarkable that this triplet 
of red galaxies lies about 1.5' (i.e., 760 kpc) eastward of the cen- 
tre of the high density region and X-ray emission described by 
Castellano et al. ( 2007). This may be an indication that the clus- 
ter is not yet relaxed. 

We also show in Fig. [2l the density of z ~ 1.6 galaxies based 
on photometric redshifts in the GOODS-MUSIC catalogue as 
filled contours (this density map is comparable to that presented 
by[Castellano et al. 2007). This density is computed by counting 
the number of neighbouring galaxies at z ~ 1.6 within a radius 
(projected on the sky) of about 300 (physical) kpc (roughly the 
scale of galaxy groups). To normalize this number and thereby 
correct for edge eff'ects and the presence of bright stars, this 
value is divided by the total number of galaxies in the cata- 
logue within this radius. Contours are subsequently drawn us- 
ing triangulation to obtain a regular grid. This process involves 
some extrapolation at the edges of the grid that produces obvi- 
ous artifacts in the density maps, which should be disregarded 
by the reader. In Fig.[2l the redshift interval of members is taken 
to be 1.43 < Zphot < 1-77, which is based on a comparison of 
photometric redshifts within the MUSIC catalogue with the red- 
shifts of spectroscopically-confirmed members. The solid and 
dashed contours are based on the photometric redshifts in the 
GMASS catalogue, using the same redshift interval as above 
but including spectroscopically-confirmed galaxies within the 
interval and excluding galaxies with spectroscopic redshifts out- 
side the interval. The result is a similar but not identical density 
map. Noteworthy is the more pronounced concentration of pho- 
tometric members, which are defined to be members according 
to photometric redshifts, in the western part of the field, which is 
close to the three massive, red z = 1.6 galaxies confirmed by the 
K20 survey ( Cimatti et al. 2004, just to the west of the GMASS 
field). The circle on this figure indicates the extent of the high- 
density region defined in Sect. 12. 1[ The radius of this circle is 
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Fig. 2. Density maps of galaxies at z = 1.6 in the GOODS-South field. Filled (solid and dashed) contours indicate the density 
of 1.43 < Zphot < 1-77 galaxies based on the GOODS-MUSIC (GMASS) catalogue photometric (and spectroscopic) redshifts. 
Contours indicate 0.0, 0.25, 0.5, and 0.75x the maximum overdensity above the median density ofz= 1.6 galaxies. This maximum 
overdensity is 1.9 (4.6) times the median density. Large dots (crosses) indicate the positions of GOODS-MUSIC (GMASS) galaxies 
in the redshift range. Small dots indicate galaxies in the GMASS catalogue outside the redshift range. CI 0332-2742 galaxies with 
confirmed redshifts are indicated by additional symbols according to their morphology: circles for elliptical, squares for spiral, and 
triangles for irregular galaxies. The elliptical galaxies have early-type spectra, except for the two indicated with double circles, 
which have intermediate-type spectra. All others have late-type spectra. Diamonds indicate galaxies with redshifts confirmed by 
ESO GOODS spectroscopy to be in the range 1.600 < Zspec < 1.622, but are not in the GMASS sample. 
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Fig. 3. The 2C long triplet of z = 1.61 early-type galaxies (as indicated by circles). Shown are a three-colour (F606W, F775W and 
F850LP) image and the individual ACS bands, provided by the publicly available HST GOODS observations. The colour scaling 
covers the range from to 25 cr, and to 40cr, respectively, where cr is the dispersion around the mean background level (which is 
zero by default here). The images are 21'' on a side. Note the much bluer foreground (and possibly background) galaxies. 



chosen to be 1 (physical) Mpc, similar to the size of low-redshift 
clusters. The area o f this circular region is similar to the area of 
a square defined by lCastellano et al.l (l2Q07l) to contain the high- 
est surface density of z = 1.61 galaxies. The high-density region 
contains 21 of the member galaxies; the same number is present 
in the area outside this region, which is five times larger and is 
therefore referred to as the low-density region. 

Selection biases in the angular distribution of galaxies with 
spectroscopic redshifts, can be caused by bright objects prevent- 
ing the detection of faint objects, for example distant galaxies, 
within a certain radial distance. Biases are also introduced by 
the design of the spectroscopic masks, which have certain ge- 
ometrical constraints such as the smallest possible distance be- 
tween objects included in diff'erent slits. To estimate the signif- 
icance of the angular inhomogeneity displayed by the members 
with redshifts determined by GMASS, we carried out a counts- 
in-cells analysis. In each cell, we compared the total number of 
objects targeted in the slits of all six GMASS masks with the 
number of z = 1.6 galaxies observed in that cell. This analysis 
was completed for a range of cell sizes, and showed that the an- 
gular distribution of the slits used for spectroscopy was rather 
homogeneous. We therefore consider any observed clustering of 
members on the sky as true angular clustering. 

3.2. MelocWy dispersion 

A fundamental property of galaxy clusters, groups, and their pro- 
genitors is the velocity dispersion of their members. In virialized 
systems, this provides a measure of the mass contained within 
the system. A widely adopted estimator of c luster velocity dis- 
persion is the biweight statistic discussed by iBeers et aP (ll99Ql 
hereinafter BFG90). For a typical dataset of 20 to 50 galaxies, 
BFG90 showed that the biweight-scale estimator is a robust esti- 
mator of distribution spread and demonstrated its insensitivity to 
outliers. We determined the median redshift and biweight statis- 
tic of the z = 1.6 spike in an iterative way. We started the itera- 
tive process using all galaxies in the GMASS field with spectro- 
scopic redshifts in the interval 1.50<z< 1.72. The median Zmed 
for these galaxies was determined, and velocities in the cluster 
reference frame were computed as follows: 

VresU = CX Az/(1 + Zmed), (1) 

where c is the speed of light. An estimate of the velocity dis- 
persion, CTspike, of the spike was provided by the biweight statis- 
tic for all Vrest,/- Subsequently, only galaxies within the velocity 
interval -3(Tspike < Vrest,/ ^ 3(Tspike were included and the com- 
putation of the median redshift and velocity dispersion was re- 
peated until no changes in the galaxies included occurred. This 



happened after two iterations; the 42 galaxies remaining inside 
the velocity interval at this stage, were defined to be members 
of the redshift spike. Their coordinates, GMASS IDs, and other 
properties are listed in Table [T] 

Our measurements of cluster redshift and velocity dispersion 
are Zmed = 1.610 and (Xspike = 440^^^ kms"\ respectively. The 
uncertainty in (Tspike was taken to be the Icr bootstrap confidence 
interval estimated from 1000 Monte Carlo simulation^. Of the 
42 members, 32 have secure redshifts determined by GMASS 
and 7 have secure redshifts determined by ESO/GOODS. The 
remaining three redshifts are less secure, but removing these ob- 
jects from the sample does not alter our final results. Expanding 
the redshift interval to a width of ±10(Tspike does not add any 
members to the 42 found. In the second inset of Fig.[TJ we over- 
lay a Gaussian function that corresponds to the computed (Tspike- 
The skewness and kurtosis of the redshifts distribution are -0.23 
and 3.6, respectively. As is evident in the histogram, there is an 
enhancement of values, relative to a Gaussian, below the dis- 
tribution mean and the distribution i s much more heavy-tailed 
than a Gaussian (Bird & Beers|[l993l) . A Gaussian function may 
not therefore be the most suitable description of the redshift dis- 
tribution. This is also shown by the more robust estimators of 
asymmetry and tail index (Bird & Beers 1993), based on order 
statistics, which have values of -0.27 and 2.6, respectively, in- 
dicating deviation from Gaussianity with significances of >20% 
and >99%. Considering only the 21 galaxies within the high den- 
sity region (dashed histogram in Fig. [T]), we obtain a velocity 
dispersion of (Thigh = 500^|qq kms"^ As this value is, within 
the errors, consistent with the value obtained above for the 42 
members, we use the value of 500 km s"^ for the overdensity's 
velocity dispersion in the remainder of this paper. 

In redshift space, the members seem to form a bimodal dis- 
tribution (Fig.[T]), with a main peak atz = 1.610 and a secondary 
peak at z = 1.602. A k-means clustering analysis indeed shows 
that, if the presence of two redshift groups is assumed, the most 
likely partitioning is into a small group of seven to ten mem- 
bers at a mean redshift z = 1.603 and a large group of the re- 
maining members with mean redshift z = 1.611. We also fitted 
two models to the redshift distribution around z = 1.61, the first 
consisting of a single Gaussian function and the second com- 
posed of the sum of two Gaussian functions. In both cases, the 
area under the curve was fixed to be the area given by the his- 
togram of the 42 member redshifts at z = 1.6, resulting in two 
and five free fitting parameters, respectively. The model com- 



^ The computation of (Tspike and its asymmetric error bars was carried 
out employing the ROSTAT software package written and provided by 
BFG90. 
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prising the summed Gaussian functions, shown in Fig.[TJ fits the 
redshift distribution more accurately: an F-test comparing these 
two models excludes the possibility that the 32% decrease in the 
value of is solely due to chance at the 98.5% level. 

The two peaks might be populated by spatially- 
separated groups. W e applied a statistical test proposed by 
iDressler & ShectmanI ([l98 8) to quantify possible spatial sub- 
structure using the positions of the galaxy members and their 
recessional velocities, but foun d no significant evid ence that 
the overdensity is substructured. iPinkney et al.l (Il996i) evaluate 
statistical tes ts for substruct ure in clusters of galaxies, inc luding 
the Lee2D (iFitchettI '1988') and Lee3D (Pinkney et al.) tests. 
They find that these t ests are less sensitive t o subclustering than 
the test introduced by lDressler & ShectmanI but have the advan- 
tage of being insensitive to non- sub structure configurations that 
appear as substructure in other tests (such as elongation and a 
velocity dispersion gradient). We applied both tests to our data 
and found marginal evidence of substructure using the Lee2D 
test and no evidence using the Lee3D test. The tests infer a 
consistent angle for the line, on which the projected member 
positions show the largest contrast at ~ 75° w.r.t. the line of 
constant declination. As this angle is almost identical to the 
position angle of the GMASS field plotted on Fig. [51 we do not 
indicate this angle separately. We note that the six masks used 
for the GMASS spectroscopy all had diff'erent position angles 
and the similarity between the two angles described above is 
considered to be a coincidence. Finally, the two -dimensional 
two-sided K-S test (iFasano & Franceschinilll987 ) also showed 
that the angular distribution of the eleven 1.600 < z < 1.608 
galaxies does not diff'er significantly from that of the 3 1 galaxies 
in the interval 1.608 < z < 1.622. The fact that no clear 
separation between the two redshift groups is seen in projection 
does not exclude the existence of substructure along the line of 
sight. 

For a relaxed cluster, the virial radius is similar to R200, the 
radius inside which the density is 200 times the critical den- 
sity. Using the redshift dependence of the critical density and 
the virial mass to relate the line-of-sight velo city dispersion to 
the cluster mass, R200 can be expressed as (see lFinn et al.ll2005l) 



^200 = 2.47 ■ 



■.h-^Mpc. (2) 



lOOOkms-i yil^ + Q^(l+^y 

In addition, the virial mass of a cluster can be expressed as 
a function of velocity dispersion using the above equation. 
Ass uming that the vi rial radius is equal to 7?2oo, one obtains (see 
also lFinnet aDl2005h 



and 



3 cr R200 



Mvir = 1.7X10^^( r\ 

VlOOOkms-i; 



loookms-i; + 11^(1+^)3 



(3) 



/i7o'Mo.(4) 



Although we do not know whether the structure under study here 
is virialized, we do apply these equations to derive an indication 
of the virial radius and mass, obtaining (for cr = 500 km s"^) 
R200 = 0.5 Mpc and M^i, = 9 xlO^^ M©. 

3.3. Galaxy overdensity and mass contained 

The galaxy overdensity in redshift space is defined to be the 
number of galaxies found in a certain redshift interval, di- 
vided by the number of galaxies expected in this interval mi- 
nus one. To determine a galaxy overdensity in redshift space. 



one has to assume an expected number of galaxies in this in- 
terval. It is obviously impossible to estimate the expected num- 
ber of galaxies without introducing any biases such as those dis- 
cussed in Sect. 12. 3[ In the absence of a superior method, we as- 
sume a flat n(z) distribution and employ the 108 galaxies with 
spectroscopically-confirmed redshifts within the GMASS sam- 
ple in the range zi = 1.400 < z < 1.900 = Z2- The lower limit 
is set by the selection criteria for GMASS spectroscopic obser- 
vations and the upper limit is based on the appearance of the 
redshift distribution shown in Fig. [H which declines at z > 2. 
Indeed the number of galaxies expected varies by less than 5% 
for 1.75 <Z2< 2.00. For the redshift interval 1.601 <z< 1.621 
with Az = 0.02, 4.32 galaxies would be expected but 40 are ob- 
served (subtracting the two galaxies that define the boundaries 
of this interval). Therefore, we derive an estimated overdensity 
of 40/4. 32- 1 = 8.3+1.5, where the uncertainty is based on 
Poisson statistics. We note that restricting the redshift interval of 
the peak to 1.601 < z < 1.615, i.e., excluding a single galaxy, 
leads to an estimated overdensity of 39/3.02 - 1 = 11.9+2.2. 
In addition, excluding the redshift interval 1.572 < z < 1.670 
from the n(z) count, which includes all member galaxies and ad- 
ditional empty redshift space, the estimated overdensity would 
be 40/3.28 - 1 = 11.2+2.2. In the following, we conservatively 
adopt an overdensity of 8.3+1.5. 

To estimate the total mass contained within the volume oc- 
cupied by the spike, we need to measure this volume and the 
matter overdensity, which is related to the galaxy overdensity by 
the bias factor 



1 + /?^mass = C(l + ^gal,obs) 



(5) 



(e.g. JSteide l et al.i ri99 8j). The bias factor is plau sibly 1.6, as ex- 
trapola ted from measurements up to z = 1.5 by iMarinoni et al.l 
(l2QQ5h. which is lower than the range of 2-4 assumed by 
ISteidel et"aD (11998 ), who originally introduced the method fol- 
lowed below. Assuming this bias factor, the matter overdensity 
then equals ^mass = 2.45 ^^3^. Here, we assume the observed 
galaxy overdensity ^gai,obs = 8.3 + 1.5, which implies a value 
of 0.53Iq q4 for the factor C and takes into account the redshift- 
space distortio ns caused by the peculiar velocities of the galaxies 
(ISteidel et al.l ll998). In the transverse direction, it is difficult to 
estimate the spatial extent of the structure because the GMASS 
field does not cover the entire overdensity. We a dopt two ap- 
proaches here: in the fi rst approach, w e follow Steidel et al.J 
(1998) and similar work d Venemans et al.' 2007; Ota et al. 2008), 
by taking the entire field sampled, while in the second, more 
conservative and maybe more accurate, approach, we take the 
localized high density region only as defined by a 1 Mpc (phys- 
ical) radius from the central galaxy. In the first case, from the 
observed redshift depth of the overdensity of Sz = 0.016, ex- 
cluding the single galaxy at z = 1.621, and taking into ac- 
count that Vtrue = Kbs/C, the mass enclosed in this volume is 
Pm V(l + ^m) = 5.9 + 0.9 X lO^'^M©. The volume considered here, 
which is equivalent to a sphere of a ~ 10 Mpc comoving radius, is 
large, corresponding to a physical size of -3.8 Mpc at z = 1.61. 
In the second case, the physical volume of a sphere has a physi- 
cal radius of 1 Mpc, but the overdensity is almost six times larger. 
The mass contained in this smaller volume is therefore a fraction 
0.018x6 of that stated above, or 6.4 + 1.0 x lO^^M©. 

The total stellar mass contained within the 21 galaxies inside 
the high-density region equals ~ 6 x 10^^ M© (the 21 galax- 
ies inside the low density region add only another ~ 2 x 10^^ 
Mo). Applying a total-to-stellar-mass ratio of 100, appropriate 
for clust ers with a velocit y dispersion in the range 425-800 
kms"^ (iBalogh et al.l l2007h , in a similar way to the approach 



8 



J. Kurk et aL: GMASS ultradeep spectroscopy of galaxies at z ~ 2 



of iMcCarthy et alJ (l2007l) for a structure at z = 1.5, we obtain a 
total mass of ~ 6 x 10^^ M©, which is consistent with the mass 
obtained above for the sphere circumscribing the high-density 
region. 

3.4. X-ray emission 

ICastellano et alJ (l2QQ7l ) reported three X-ray sources within the 
high-density region, one coinciding with the peak of their galaxy 
density contours and including a late-type galaxy of a peculiar 
morphology. We identify this galaxy with GMASS 2540, which 
indeed has a redshift z = 1.613, measured using GMASS spec- 
troscopy. In addition, several other galaxies with spectroscopic 
redshifts at z = 1.6 from GMASS ar e located within the soft 
X-ray contours displayed in Fig. 2 of ICastellano et al.L namely 
GMASS 2180, 2286, and 2251 at z = 1.608, 1.604, and 1.609, 
respectively. 

The CD FS was observed by Chandra for 1 Msec 
(iRosati et al.l l2002 b). which allows detection of extended 
emission down to 1 x 10"^^ ergcm"^s"^ in the 0.5-2 keV 
band for relatively compact sources. More recently, t hese 
observations were extended to 2 Msec (iLuo et al.l l2008b . In 
the GMASS region, the Chandra PSF FWHM is 2-3 arcsec 
and the median effective exposure time at the positions of 
the overdensity members is 1.56 Msec. For sources with an 
extent of 5-10 arcsec, which are easily resolved by Chandra, 
the 2 Msec data have a detection limit (S/N ^ 3) of 6 x 10"^^ 
erg cm"^ s"^ in the soft band. Given the size of the GMASS field 
and the cumulative X-ray cluster number counts presented in 
iRosati et aP (l2002ah , at least ~ 4 low X-ray luminosity clusters 
(at any redshift) are expected to be detected. The presence of a 
galaxy group or cluster in this field is therefore quite likely. 

At z = 1.6, this X-ray sensitivity implies a limit to the 
X-ray luminosity of 1 x 10^^ erg s"^ in the observed frame 
and 0.5-2 keV band. For a T= 2 keV group, the K-correction 
amounts to 9 and the bolometric correction to 2.3. Therefore, 
the Cha ndra observations i mply a limit of Lx,boi - 2.1 x 10^^ 
erg s~^ iMahdavi & GeUeJ (l200l]) provide an empirical relation 

= 5 X 10^^ erg s"^ ((t/330)'^ '^, which has significant scatter, 
in particular for groups. For the measured velocity dispersion of 
the redshift spike of 500+100 km s"^ one would expect an X- 
ray luminosity 3^2 10"^^ erg s"^ ; this is consistent with our limit 
because of both the large scatter in the Lx - cr relation, and the 
error in the velocity dispersion, which may have been overesti- 
mated if the overdensity is not yet virialized. An extended X-ray 
source is detected neither in the vicinity of any ETG in the north 
of the field, nor around the triplet clump. In addition, by stacking 
X-ray photons in the 0.5-2 keV band around the seven ETGs, for 
a total eff'ective exposure time of 1 1 .3 Msec, no signal is detected 
down to 0.7 X 10"^^ erg cm"^ s"^ 

3.5. Colour-magnitude diagram 

Although more pronounced in galaxy clusters, the field popula- 
tion of galaxies up to at least z ~ 1 shows a bimodal distribution 
in colour, with red galaxies on the passive side and blue galax- 
ies on the active side, separate d by a green vall ey of galaxies in 
transition from blue to red (e.g. jBell etal .120041) . The emergence 
of colour bimodality at z ~ 2 is the subject of an other paper em- 
ploying the GMASS spectra ("Cassata et aLll2008h . Within cluster 
populations, the bimodality is detected as a concentration of red 
(early-type) galaxies on a line in a colour-magnitude diagram 
using colours that straddle the 4000 A break, which is referred 



to as the red sequence (RS). This sequence is more pronounced 
in regions of high galaxy density (i.e., clusters and groups) be- 
cause galaxies appear to evolve more rapidly in such an environ- 
ment; Thoma s et al.l (l2005l) . for example, found that most of the 
star-formation activity in early-type galaxies probably occurred 
at 3 < z < 5 in high, and at 1 < z < 2 in low density environ- 
ments. The RS was found to be present in a cluster at z = 0.83 
( Tran et al. 2007) , has been used to find clusters up to z ~ 1.4 
( Gla dders & Yed l2000), and there is evidence that the bright end 
(M:, > 10^^ Mo) of the RS is populated up to z ~ 2 in massive 
cluster progenitors (iKodama et al .112007b . 

For z = 1 .6, the 4000 A break is shifted to 10 400 A, between 
the z and / bands. We plot a z-/ versus / diagram in Fig.lH indi- 
cating the member galaxies by large symbols (defined according 
to spectral class: circles for early types, diamonds for late types, 
and both for intermediate types). Also indicated, with small cir- 
cles, are objects without a spectroscopic redshift but with a pho- 
tometric redshift in the range 1.50 < Zphot < 1-70. There appear 
to be two populations of spectroscopic members: a blue popula- 
tion around z - / = 0.6, and a red population with z - / ~ 1.5. 
Clearly the early and intermediate types (derived from the spec- 
tra) have the reddest z - J colours: all are at z - / > 1.3. The 
figure shows that, although there are many member candidates 
present at z - / > 1.3, there is also some contamination from 
objects (as indicated by the crosses) with spectroscopic redshifts 
outside the range 1.600 < Zspec < 1.622. 

4. Photometric, morphological and spectral, galaxy 
properties in spike and field 

In this section, we compare the galaxy properties of samples of 
member and field galaxies. The properties include the galaxy 
morphology, spectral type, fitted SE P, and severa l properties 
derived froni this SED. We refer to Cassata et al. (2008) and 
ICimatti et al.l (2008) for more details about the morphological 
analysis and SED fitting. 

4.1. Morphological analysis 

The morphological analysis was performed independently by 
two of us (P C and GR) in the ACS F850LP-band image 
(Cassata et al. 2008). Three classes were recognized by eye: 
(1) elliptical, (2) spiral, and (3) irregular. For the faint galaxies 
(close to z = 1.6), the two classifiers agreed in more than 70% 
of the cases. The discrepancies mainly occur for galaxies close 
to the boundary between spirals and irregulars, whereas agree- 
ment was obtained for 90% of the early -types. The final cata- 
logue has been produced by a reconciliation process between 
the two classifiers for the remaining 30% of objects. We also 
determined morphological parameters by fitt ing the galaxy light 
profiles (for details of the fitting process, see lCassata et al.h . We 
compared the classification by eye with the machine-determined 
parameters, obtaining the following results: 90% of the early- 
type galaxies (12/13) have a Sersic index > 2, while 88% of late- 
types (63/72) have a Sersic index < 2. 

The morphological classification was carried out using the 
ACS F850LP-band image, i.e., using the rest- frame UV light 
of galaxies at z = 1.6. Therefore, surface brightness dimming 
and morphological k-correction may be aff'ecting the morpho- 
logical classification. On the one hand, ellipticals are expected 
to remain symmetrical in the UV, but fainter because of the k- 
correction eff'ect. On the other hand, spirals can appear morpho- 
logically later because knots of star formation brighten in the UV 
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Table 1. Spike galaxy redshifts, photometry, properties derived from fitted SEDs, and morphological classifications 



ID 



R.A. & Dec. 



4.5yt/m 




MS 



365 
648 
781 
1254 
1282 
1295 
1380 
1495 
1556 
1667 
1691 
1708 
1808 
1871 
1979 
1988 
2055 
2081 
2111 
2113 
2142 
2148 
2180 
2196 
2251 
2286 
2341 
2352 
2355 
2361 
2368 
2438 
2454 
2493 
2540 
2543 
2550 
2603 
505 
1081 
1399 
1413 



03 32 27.80 
03 32 29.81 
03 32 17.71 
03 32 20.17 
03 32 33.86 
03 32 28.59 
03 32 25.25 ■ 
03 32 35.36 
03 32 37.74 
03 32 40.99 
03 32 31.90 
03 32 23.12 
03 32 26.15 
03 32 26.82 
03 32 24.64 
03 32 24.02 
03 32 26.77 ■ 
03 32 29.86 
03 32 27.94 
03 32 22.00 
03 32 23.54 
03 32 36.30 
03 32 29.56 
03 32 36.67 
03 32 29.48 
03 32 29.99 
03 32 17.52 
03 32 33.88 
03 32 14.32 
03 32 26.05 
03 32 17.10 
03 32 26.41 
03 32 28.91 
03 32 38.51 
03 32 30.33 
03 32 35.92 
03 32 30.08 
03 32 27.85 
03 32 43.37 
03 32 17.51 
03 32 41.66 
03 32 24.66 



27 48 12.0 
■27 47 19.7 
27 47 02.9 
27 45 49.3 
27 45 42.6 
■27 45 42.2 
27 45 29.0 
27 45 12.6 
■27 45 05.5 
27 44 50.2 
27 44 45.0 
27 44 42.2 
■27 44 33.3 
.21 44 24.7 
27 44 07.8 
27 44 08.3 
27 43 58.1 
27 43 54.8 
■27 42 45.7 
27 42 43.5 
27 42 49.3 
27 42 49.5 
27 42 56.0 
27 42 58.5 
■27 43 22.0 
■27 43 22.6 
27 43 36.6 
27 42 04.1 
■27 43 32.9 
27 42 36.6 
27 43 41.9 
■27 42 28.4 
27 43 03.6 
■27 42 28.0 
27 42 40.3 
27 42 41.0 
■27 42 12.2 
27 43 05.7 
■27 47 43.5 
27 46 14.0 
■27 45 25.6 
27 45 21.6 



1.609 
1.609 
1.605 
1.610 
1.621 
1.611 
1.612 
1.611 
1.607: 
1.613 
1.612 
1.607 
1.609 
1.613 
1.611 
1.610 
1.611 
1.601 
1.610 
1.613 
1.610 
1.609 
1.608 
1.614 
1.609 
1.604 
1.602 
1.604 
1.610 
1.609 
1.612 
1.615 
1.602 
1.607 
1.613 
1.612 
1.601 
1.612 
1.610: 
1.610 
1.615 
1.614: 



2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 
1.0 
2.0 
2.0 
1.0 
2.0 
1.0 
2.0 
1.0 
2.0 
2.0 
1.0 
1.5 
2.0 
2.0 
2.0 
2.0 
2.0 
1.5 
2.0 
2.0 
2.0 
2.0 
2.0 
2.0 



23.10+0.02 
22.20+0.01 
22.54+0.02 
22.21+0.02 
21.50+0.01 
23.08+0.03 
22.30+0.02 
21.20+0.01 
21.03+0.01 
22.51+0.02 
21.11+0.01 
20.81+0.01 
22.61+0.02 
22.66+0.02 
20.95+0.01 
22.77+0.03 
22.21+0.02 
22.09+0.01 
20.64+0.01 
20.37+0.01 
22.22+0.02 
19.80+0.01 
21.00+0.01 
20.34+0.01 
20.64+0.01 
20.82+0.01 
22.47+0.02 
20.09+0.01 
20.92+0.01 
20.40+0.01 
21.15+0.01 
20.87+0.01 
20.80+0.01 
20.88+0.01 
21.40+0.01 
20.26+0.01 
22.03+0.01 
21.74+0.01 
22.86+0.02 
21.03+0.01 
22.71+0.02 
21.50+0.01 



24.4+0.1 
23.9+0.1 
23.5+0.1 
23.6+0.1 
23.9+0.1 
24.5+0.1 
24.1+0.1 
24.0+0.1 
24.8+0.2 
25.1+0.2 
23.6+0.1 
25.0+0.2 
23.9+0.1 
23.4+0.1 
23.2+0.1 
24.0+0.1 
24.8+0.2 
24.1+0.1 
24.1+0.1 
24.4+0.1 
24.0+0.1 
23.4+0.1 
23.4+0.1 
23.8+0.1 
23.7+0.1 
24.3+0.1 
24.0+0.1 
23.4+0.1 
24.2+0.1 
23.9+0.1 
24.0+0.1 
23.5+0.1 
24.0+0.1 
23.3+0.1 
23.7+0.1 
24.6+0.1 
23.9+0.1 
23.6+0.1 
24.5+0.1 
23.4+0.1 
24.6+0.1 
23.8+0.1 



23.9+0.1 
23.2+0.1 
23.3+0.1 
23.4+0.1 
23.1+0.1 
23.7+0.1 
23.4+0.1 
23.3+0.1 
23.6+0.1 
24.5+0.1 
22.9+0.1 
23.8+0.1 
23.5+0.1 
23.0+0.1 
22.7+0.1 
23.6+0.1 
23.9+0.1 
23.4+0.1 
22.4+0.1 
23.2+0.1 
23.3+0.1 
21.8+0.1 
22.7+0.1 
22.2+0.1 
22.6+0.1 
22.9+0.1 
23.4+0.1 
22.1+0.1 
22.7+0.1 
22.3+0.1 
23.3+0.1 
22.9+0.1 
23.4+0.1 
22.6+0.1 
22.9+0.1 
22.8+0.1 
23.3+0.1 
23.1+0.1 
23.9+0.1 
22.8+0.1 
24.1+0.1 
23.2+0.1 



Spectroscopic redshift, colon indicates uncertain redshift 
Origin of redshift: (1) GMASS (2) ESO/GOODS 

Spectral class: (1.0) early type, (1.5) early type with signs of star formation, (2.0) late type 
Logarithm of stellar mass in solar masses 
SFR in solar mass per year 
Age in Gyr 

Morphology: (1) Elliptical (2) Spiral (3) Irregular 



and disks become fainter because of the surface brightness dim- 
ming. The combination of these two eff'ects should not aff'ect our 
visual classification, at least in separating early- and late-types. 
As the galaxies in our samples are observed in almost the same 
rest-frame UV light (the central rest-frame wavelength diff'ers 
from 394 to 338 nm for redshifts between z = lA and z = 1.8, 
respectively, while the filter has a width of 47 nm in the rest- 
frame), these phenomena do not introduce a diff'erential eff'ect 
on the classification of galaxies at z = 1.4 to z = 1.8. 



4.2. SED fits 

The properties of the stellar populations in each galaxy are de- 
termined from the SED fitting analysis for the spectroscopic 
redshift determined. The results of the SED fitting are, there- 
fore, far more secure than for one where the redshift is unknown 
and therefore a free parameter in the fitting process, as is the 
case for the determination of pho tometric redshift s. We adopt 
the libraries of synthetic spectra by'Maraston' ("2005', M05), who 
employs the Kroupa (2001) initial mass function (IME) and in- 
cludes the thermally pulsing asymptotic giant branch (TP-AGB) 
phase of stellar evolution. This phase is the dominant source of 
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Fig. 5. Comparison of galaxy properties for the galaxies in the redshift range 1.4 < z < 1.8. From left to right and top to bottom: 
SFR, mass, sSFR, and age. Symbols are for galaxies within the high density region (filled circles), low density region (open circles), 
and outside the peak (boxes). Spectroscopic early-type galaxies are indicated by stars. The dashed lines in the fourth panel indicates 
the age of a galaxy formed with a short burst of star formation at z = 2.0 and at z = 2.8. 



bolometric and near-IR radiation for a simple stellar population 
in the age range from 0.2 to 2 Gyr 

We adopt exponentially declining star-formation histories, 
i.e., SFR = (M/T)exp(-t/T) with r =0.1, 0.3, 1, 2, 3, 5, 10, 
15, and 30 Gyr, and, in addition, the case of a constant star- 
formation rate. Extinction is treated as a free parameter in the 
optimization, by adopting the extinction curve of Calzetti et al. 
(2000). We assumed solar metallicity for all the models. The fit- 
ting procedure selects the template spectrum that minimizes the 



X^, and therefore provides a value for each of stellar population 
properties: the age of the best-fit model, the ^-folding time of the 
SFR T, the extinction Ay, and the stellar mass. In addition, the in- 
stantaneous SFR is computed from the age and r. During the fit- 
ting, only the observed bands corresponding up to the rest- frame 
^^-band (Arest < 2.5yum) were used to avoid any dust-emission 
contamination. While the details of t he procedure and its un- 
certainties are extensively discussed in lPozzetti et al.l (l2007l) . we 
recall a few points here. The median formal statistical uncer- 
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Fig. 4. Colour-magnitude diagram of all objects in the GMASS 
catalogue (excluding those undetected in one of the two bands) 
showing z - J colour vs. / magnitude. Large symbols indicate 
galaxies at 1.600 < Zspec < 1.622, small circles indicate galaxies 
at 1.50 < Zphot < 1-70 (not including those with Zspec < 1.600 
or ^spec > 1.622). Galaxies with spectroscopic redshifts out- 
side the range 1.600 

^ ^spec ^ 1.622 are indicated by crosses. 
Spectroscopic early types are indicated by circles (5), interme- 
diate types by a both circles and boxes (2), and late types by 
boxes only (35). The solid (red) line indicates a theoretical red 
sequence observed at z = 1 .60 for elliptical galaxies with a range 
of masses (as indicated in logarithmic solar masses) formed at 
Zf = 3.0. At the bottom of the plot, typical errors in colour are 
indicated for each / magnitude. In addition, the objects to the 
right of the diagonal dashed line have data of signal-to-noise ra- 
tios below three. To the right, a histogram is displayed for all 
member galaxies (open histogram) and for members within the 
range 21.5 </< 23.0 (indicated by dashed vertical lines in the 
colour-magnitude diagram) where the confirmed red sequence 
galaxies are located (filled histogram). 



tainties, derived from the width of the probability distributions 
in the fitting procedure, for each stellar population parameter 
(mass, age, SFR, Ay, r) are appr oximately 10-20%. U sing ex- 
tensive Monte Carlo simulations, Pozzetti et al.l (l2007h demon- 
strated that the overall internal accuracy of the measured stellar 
masses is ~ 0.2 dex. In addition to this typical internal error, 
systematic errors should also be mentioned. One source of sys- 
tematic error is the addition of secondary bursts to a continuous 
star formation history, which produces systematically higher (up 
to 40% on average) stellar masses. In addition, we note that the 
synthetic spectra including the TP-AGB phase (lMarastonll2005[ 
Chariot & Bruzual, in prep.) produce, instead, a shift of ~ 0.1 
dex towards low er stellar masses than models without this im- 
portant phase (Bruzual & Chariot 2003), at least for our sample. 
Finally, the uncertainty in the absolute value of the stellar mass, 
owing to assumptions about the IMF, is within a factor of two 
for the IMFs typically adopted in the literature. 



The galaxies with log(sSFR) = -7.9 are aff'ected by the limit 
set to be the minimum age of 9 x 10^ yr. Since an exponential 
star formation history is assumed, the sSFR estimate is related 
to the age in the following way: sSFR oc exp -t/r. 

There are no strong correlations between the derived prop- 
erties, apart from the following classical ones: (1) the age- 
metallicity correlation, and (2) the age-reddening correlation for 
red galaxies. The first relation is not applicable here because we 
do not derive the metallicity (it is fixed to be solar). The second 
one may be a concern as it can cause a degeneracy between dust- 
free galaxies with old stellar populations and dusty galaxies with 
young populations. However, SED fitting is able to break this de- 
generacy if sufficient bands are included around the (r edshifted) 
4000A break (see for example iPozzetti & Mannuccil l^OO. for 
EROs, using the RJK bands). In the case ofz= 1.6 galaxies, the 
4000A break is shifted to 1 fim. As our SED fitting of galaxies at 
z = 1.6 includes bands up to 6.5 yum, this degeneracy should not 
be a major concern. 



4.3. Comparison of field and overdensity samples 

To compare the properties of the member galaxies with those of 
field galaxies, we selected a sample of 43 galaxies outside the 
peak, that is 24 in the redshift interval 1.416 < z < 1.598 and 
19 in the interval 1.624 < z < 1.840, which we refer to as the 
field sample. The field galaxies were selected from the GMASS 
sample and have therefore similar information available, derived 
in an identical way, to that for the spike galaxies. Of the 43 field 
galaxies, 25 have redshifts determined by GMASS only, and 8 
have uncertain redshifts. We further subdivide the 42 galaxies 
within the redshift peak into those within the (spatially) high- 
density region defined in Sect. [T] (21 galaxies), and those outside 
(also 21, low-density region). 

For an unbiased comparison of derived galaxy properties, it 
is necessary to know whether there are any selection eff'ects re- 
lated to the galaxy redshifts. The mass limit is mainly defined by 
the m4.5^m < 23.0 limit, which sets a completeness limit of about 

= 3 X lO^Mo at z < 1.8 (for the M05 templates used). There 
are six galaxies with masses below this limit within and six out- 
side the redshift spike. We set the fitted mass of these galaxies 
to log(M*) = 9.5Mo when we compare the sample masses. All 
SFRs below log(Mo yr"^) = -0.25 are set to this lower limit be- 
cause they are not detected in the U' band and SFRs below this 
value can therefore not be measured accurately. A lower limit 
to the age (of 9 x 10^ yr) has been set by the fitting program. 
This enforces an upper limit to the sSFR because the stellar mass 
is computed by assuming an exponential star formation history 
over the lifetime of the galaxy. In Sect. l2.3[ we demonstrated that 
the [O II] detectability is not responsible for the high density of 
spectroscopically confirmed galaxies at z = 1.6. The contribu- 
tion of [Oil] detected galaxies to the sample at z > 1.6 is nev- 
ertheless much smaller, which might cause a bias in the sense 
that fewer galaxies with high star-formation rate are selected. 
However, the upper-left plot in Fig. [5] indicates that this is not 
the case. In contrast, the average SFR (derived from SED fitting) 
of spectroscopically confirmed z > 1.6 galaxies is higher than in 
the spike. We do not expect a redshift bias to be present for any 
of the other properties measured. 

We compare the following properties: stellar mass, star for- 
mation rate (SFR), specific SFR (sSFR), i.e., SFR/mass, age, rest 
frame B- 1 colour, spectroscopic galaxy class, and morphologi- 
cal class (see Table [TJ. 
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Table 2. Mean and error in the mean of galaxy properties Morphological (left) and spectroscopicai (right) classification 




High Low Spike Field 

Prop Unit Mean (Standard deviation) 

M^s log Mo 10.2 (0.3) 9.8 (0.2) 

SFR log Mo yr-i 1.0 (0.2) 1.4 (0.1) 

sSFR logyr-i -9.2(0.2) -8.4(0.1) 

Age lO'^yr 0.6(0.1) 0.2(0.1) 

B-I 0.9(0.1) 0.7(0.1) 



Table 3. Results of the K-S tests comparing galaxy properties 



Samples^ 


High-Low 


High-Field 


Low-Field 


S-Field 


Prop 


K-S 


K-S 


K-S 


K-S 


Mass 


0.029 


0.092 


0.462 


0.843 


SFR 


0.304 


0.098 


0.932 


0.506 


sSFR 


0.004 


0.010 


0.936 


0.204 


Age 


0.029 


0.018 


0.688 


0.473 


B-I 


0.071 


0.462 


0.418 


0.871 


Class 


0.583 


0.023 


0.106 


0.741 


Morph 


0.071 


0.011 


0.523 


0.923 



^ High and low density region, in spike (S) and field (F) 



The spectroscopic class is determined by the presence of ab- 
sorption lines only (1.0), emission lines only (2.0), or absorption 
lines with some emission lines present (1.5), in the GMASS and 
ESO/GOODS spectra. This classification corresponds roughly to 
passive, early-type (1.0), actively star-forming, late-type (2.0), 
and intermediate-type (1.5) galaxies. We note that the spectral 
type correlates well with the observed morphology, that is early- 
types have compact bulge-like morphologies, while late-types 
often show evidence for disk-like or irregular morphologies (see 
Figs. lA.ll to IA.4b . There is also very good agreement between 
the properties deriv ed from the SED fitting and the spectra. 
ICimatti et al.l (l2008h show that the galaxies in a sample of high- 
redshift (supposedly) early-type galaxies from GMASS are red, 
have spheroidal morphologies and spectra consistent with a pop- 
ulation of old, nearly passively evolving, stars. In addition, the 
photometric SED fitting of these galaxies gives consistent re- 
sults. 

Some properties of the spike and field samples are shown 
in Fig. \5\ while the results of K-S tests, which determine the 
probability that the properties of the various samples are drawn 
from the same distribution, are listed in Table [3] The values in 
this table are probabilities, that is a lower number indicates a 
more significant diflTerence in values for the various samples. The 
samples compared are high versus low density, high versus field, 
low versus field, and spike versus field, respectively. The mean 
and the error in the mean of these properties are listed for each 
sample in Table [2l 

There are diflTerences for all properties listed between the 
high-density sample and the low-density or field sample, al- 
though some are only marginally significant. We indicate the 
probabilities of below 3% in boldface in Table [3l In partic- 
ular, the diflTerences in sSFR and morphology are significant. 
These are probably driven by the presence of passively evolving 
galaxies with early type morphology in the high-density region, 
which are not present in the low-density region. 

The mean diflTerences between the high-density region and 
the field amount to 0.2, 0.5, and 0.7 dex in mass, SFR, and 
sSFR, respectively, and 0.3 Gyr in age. Although there are sev- 
eral galaxies with masses M > 10 ^^-^M© in the field sample. 




Fig. 6. A comparison of the galaxy morphologies and spectral 
classes between the three samples (high and low density region, 
and field). 



there are twice as many galaxies with M > 10 ^^-^M© in the high- 
density sample as in the field sample, which is twice as large (six 
against three). The massive member galaxies are also those de- 
termined from their spectra to be early types and have the lowest 
sSFRs. Within the spike, there are seven galaxies of early and 
intermediate spectral class, while outside there is only one such 
galaxy. We note that the more massive galaxies have higher ages. 

There is a notable diflference in the angular distribution of 
morphologically-classified elliptical galaxies in the spike: of the 
eight elliptical galaxies, six are inside the high-density region 
(see also Fig. [2]). Figure [6] shows a comparison between the three 
samples of the morphologies and spectral classes of the galaxies. 

We also studied the correlation between the z-J colour, spec- 
troscopic class, and morphology and the smaller-scale galaxy 
density, as computed for and shown in Fig. O for each galaxy. 
However, we did not find convincing evidence of any significant 
correlations. 

Another interesting comparison is listed in Table (H where 
values below 0.001 are printed in boldface. We display the prob- 
abilities that the following samples have properties derived from 
the same distribution: blue-member versus red-member galax- 
ies, blue-member versus blue-field galaxies, red-member versus 
red-field galaxies, and blue-field versus red field galaxies, re- 
spectively. As expected, there are significant diflTerences between 
most properties (mass, sSFR, age, and morphology) of the blue 
and red members. The properties of blue galaxies inside and out- 
side the spike do not diflTer significantly, but the red members in- 
side the spike do have lower SFRs and sSFRs than the red field 
galaxies (the mean sSFR having an order of magnitude diflTer- 
ence). In the field, only the diflTerences in mass and spectroscopic 
class are significant between the blue and red galaxies. This sug- 
gests that the red galaxies inside the overdensity have evolved 
more rapidly than those outside, while the blue galaxies appear 
to have had a similar evolution inside and outside the overden- 
sity. 



10.0 (0.2) 10.0 (0.2) 

1.2(0.1) 1.5(0.1) 

-8.8(0.2) -8.5(0.1) 

0.4(0.1) 0.3(0.1) 

0.8(0.1) 0.8(0.1) 
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Table 4. Blue and red galaxy property comparison 



Samples^ 


BS-RS 


BS-BF 


RS-RF 


BF-RF 


Prop 


K-S 


K-S 


K-S 


K-S 


Mass 


<0.0001 


0.818 


0.525 


0.0003 


SFR 


0.0031 


0.830 


0.142 


0.127 


sSFR 


<0.0001 


0.592 


0.103 


0.451 


Age 


0.0002 


0.741 


0.051 


0.598 


Class 


0.0124 


0.983 


0.103 


0.122 


Morph 


0.0006 


0.980 


0.421 


0.200 



Blue (B) and red (R) galaxies, in spike (S) and field (F) 



I 
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Fig. 7. Shown here are rest-frame composite spectra of several 
groups of CI 0332-2742 galaxies with GMASS spectra. The top 
panel shows a comparison of the composite spectra of galaxies 
in the high-density region (black) and in the low-density (grey). 
The middle panel shows a comparison of red (black) and blue 
(grey) CI 0332-2742 spectra, while the bottom panel shows the 
composite spectra of early (black) and late-type (grey) galaxies 
in CI 0332-2742. Absorption and emission features (possibly) 
identified are indicated by vertical lines below and above the 
spectra. In the middle and lower panel, the spectra are oflTset in 
vertical direction by an arbitrary amount to increase their visi- 
bility. 



4.4. Composite spectra 

Finally, in Fig. [71 we compare the rest-frame composite spectra 
of the 18 CI 0332-2742 galaxies that have GMASS spectra in 
the high-density region, with the composite spectrum of the 14 
galaxies in the low-density region. Also shown in this figure are 
composite spectra of the 1 1 red and 21 blue galaxies and of the 7 
early- and 25 late-type galaxies in CI 0332-2742. The compos- 
ites have been normalized in the interval 2300 - 2700 A, apart 
from the red and early-type galaxy composites, which were nor- 
malized in the interval 3000 - 3500 A. It is clearly noticeable that 
the composite spectrum of the galaxies in the high-density region 
is far redder and its [Oil] line equivalent width far lower (10 A 



in contrast to 40 A) than the composite spectrum of the galaxies 
in the low-density region. The same is true for the red and blue, 
and early- and late-type galaxy composites (13 A compared to 
28 A, and 2 A compared to 36 A, respectively). 

Apart from the 4000 A and 3600 A B aimer break (which are 
at /I > 1 yum for z > lA galaxies), another strong age-dependent 
feature in the UV is a break in the region 2640 - 2850 A, which is 
produced by the combi nation of sev e ral str ong absorption lines, 
including Mg 11/12800. iDaddi et al.l (l2005h named this feature 
Mguv and defined it to be the ratio of (two times) the flux at 
2625 < A< 2725 A t o the flux in t wo 100 A broad regions above 
and below this range. IPaddi et al.l show that only stellar popula- 
tions older than 0.1 (0.4) Gyr have Mguv significantly above 1.1 
(1.4). Younger stellar populations and populations with constant 
star formation (but reddening E(B-V)=1.2) have Mguv near 1.0. 
This feature is, therefore, a key signature of passively evolv- 
ing stellar populations because it is not present in young dust- 
reddened star-forming galaxies: this therefore allows us to break 
the age/dust degeneracy aflfecting red sources. The measurement 
of the Mguv index for our composite spectra clearly diflTer: the 
high-density-region, red, galaxy composite has an Mguv mea- 
surement of 1.2 in contrast to 1.1 for low-density-region, blue 
galaxies, while the early- and late-type composites have Mguv 
values of 1 .4 and 1.1, respectively. 

5. Discussion 

We have presented evidence of a structure of galaxies at z = 1.6, 
which may be a cluster in formation. The constituent galaxies 
exhibit a galaxy overdensity of at least a factor of eight and a red- 
sequence of passive galaxies, and they appear to be both more 
massive and evolved than galaxies in the field at similar redshift. 

We have compared this structure with clusters at z > 1.0, 
which have similar properties, and often less reliably measured 
spectroscopic redshifts available. The diflTerences in galaxy prop- 
erties found by ourselves between the high density region, the 
low density region, and the field are consistent with the density 
relations found at z ~ 1 by Smith et al. ( 2005) and lPostman et al.l 
(200 51), at z = 1.24 by Demar co et al.1 (l2007), and at z = 1.5 
by Fassbender et all (12008): red early-types dominate the clus- 
ter core, while blue late-types are found in the outskirts of the 
cluster. 

5.1. Velocity dispersion 

The velocity dispersion determined for the galaxies in the high- 
density region and for all galaxies in CI 0332-2742 are consis- 
tent with a val ue of 500 km s ~^. Th is va lue is about the aver- 
age va lue that iHalliday et al.l (|2004|) and iMilvang- Jensen et aP 
(12008 ) found for a sample of 26 clusters at 0.40 < z < 0.96 se- 
lected from the ESO Distant Cluster Survey (EDisCS), but lower 
than the clusters a t 0.8 < z < 1. 3 observed by the ACS/HST in 
GTO time (Postman et al.ll2005l) , most of which were first de- 
tected by means of X-ray emission. The spike at z = 1 .6 is there- 
fore, as for the EDisCS clusters, more likely to be a proge nitor 
of typical (less rich) low-redshift clusters ( Mil vang- Jensen e t al.l 
2008). 

The distribution of members of CI 0332-2742 in redshift 
space is bimodal (Fig.[T]), w ith a main peak at z = 1.610 and a 
secondary peak at z = 1 .602. lDemarco et al" I (l2007h found a sim- 
ilar bimodal redshift distribution among the 38 confirmed mem- 
bers of the cluster RDCS J1252.9-2927 at z = 1.237. In contrast 
to the latter cluster, we do not find evidence that the substructure 
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in redshift space consists of two spatially separated groups, but it 
is probable that also in the case of CI 0332-2742 the bimodality 
shows that the structure is not yet virialized. 

5.2. Red sequence 

In the colour-magnitude diagram of the observed field (Fig. |4]), 
we overplot a line that represents a theoretical red sequence at 
z = 1.60, which corresponds to the predicted magnitudes in z- 
band and /-band of the stellar populations of elliptical galax- 
ies that formed in a short time (0.5 Gyr) at Zf = 3.0 (provided 
by T. Kodam a, see Kodama et al. 1998, based on the code by 
iKodama & Arimoto 1997) . We indicate the (logarithm of) stel- 
lar masses corresponding to the positions on the red sequence. 
The early-type galaxy members are straddled about this red se- 
quence, with a scatter in z - / of 0.147^q q^3 magnitude. This 
scatter also includes the scatter produced by photometric er- 
rors in the measured colour. Subtracting this error (averaged in 
quadrature over the seven galaxies) in quadrature, we obtain an 
intrinsic scatter of 0.135 magnitudes. The tightness in the ob- 
served red sequences at low redshift is due to the large time 
diff'erence between the epoch of formation and the moment of 
observation. Observing red sequences at higher redshift, we ap- 
proach this epoch of formation and expect the scatter to increase 
due to small diff'erences in the formatio n time and evolution of 
the early-type galaxies. [Stanford et al.l C1998) found an intrin- 
sic scatter of about .05 for clusters from z = 0.1 to z = 0.8, 
and lEUis et al. 1 (1200^) measured 0.2+0.2 in a massive cluster at 
z = 0.89. The scatter that we measure at z = 1.61 is consistent 
with that expected from observations at lower redshift, taking 
into account the younger average age of the stellar populations 
at t his redshift. 

iDe Lucia et al.l (l2007l) found a clear decrease in the ratio 
of luminous (L > 0.4L*) to faint red galaxies in clusters from 
z ~ 0.8 to z ~ 0.4. This increasing fraction of faint red galax- 
ies towards lower redshift could be explained if the RS of high 
redshift clusters does not contain all of the progenitors of nearby 
RS cluster galaxies. Instead, a significant fraction of these must 
have moved on to the RS below z ~ 0.8. Although Fig. [4] appears 
to support this claim by the lack of galaxies with Zphot = 1.6 on 
the red sequence at / > 23.0, we note that, at these fainter mag- 
nitudes, we have larger errors in the photometry and photometric 
redshifts. We are unable to confirm this decrement in the num- 
ber of faint galaxies on the RS, spectroscopically, because of the 
limited completeness of our spectroscopy. 

The red sequence appears to be ubiquitous in the field and in 
clusters at low redsh ift, and is also detected at redshifts z > 1 . 
[McC arthy et al. (2007) detected one of the most distant red se- 
quences discovered so far, at z = 1.5, including at least ten galax- 
ies with very red optical-to-NIR and optical-to-MIR colours 
within a circle of 128 kpc radius, one of which has a spec- 
troscopic redshift of z = 1.51. The total stellar mass of these 
galaxies is ~ 8 x 10^^ M©, which is similar to the total stellar 
mass computed for CI 0332-2742. Six of the galaxies lie on a 
red sequence in a / - ^ CM-diagram, which is consistent with 
the colour-magnitude relation for the Coma cluster evolved to 
z = 1.5, assuming a passively-evolving SSP with Z/ = 5 (ac- 
cording to lMcCarthv et al.ll20Q7l) . 

The bimodali ty in g alaxy colours up to z ~ 2 was also found 
bvlCassata et al.1 (l2008h in the GMASS field and lGiallongo et al.1 
0005^ in the HDF. The latter also showed that the blue and red 
galaxy populations most probably had diff'erent evolutionary his- 
tories. Using a much larger sample, including -22000 galaxies 
selected over an area of 0.6 deg^ from the Early Data Release of 



the UKIDSS Ultra Deep Survey, ICirasuolo et"aD (l2007l) found 
that the colour bimodality disappears at z > 1.5. Recently, 
Williams et al. ( 2009) showed, with a large ^-selected sample 
of more than 30000 galaxies, that the bimodal distribution of 
star-forming and quiescent galaxies is still seen in a subsample 
of galaxies in the redshift range 1.5 < z < 2, but not above 
z = 2. The clear bimodality demonstrated for CI 0332-2742 
may indicate that the colour bimodality develops earlier in high- 
density regions tha n in a field samp le, such as that employed by 
ICirasuolo et al.l and I Williams et al.l . 

The relationship bet ween galaxy colour a nd environment is 
studied in more detail by Coope r et al.l (|2007, CNC07) employ- 
ing a sample of 19 464 galaxies selected from the DEEP2 Galaxy 
Redshift Survey. They measured the galaxy surface overdensity 
within redshifts bins up to z ~ 1.4 by calculating the distance to 
the third closest neighbour, and the red fraction by considering 
(U - 5)-restframe colour, which is slightly dependent on lumi- 
nosity. The overdensities sampled range up to about ten, which, 
as CNC07 state, does not cover the regime of massive clusters. 
CNC07 stated that the red fraction at z > 1.4 is not higher in 
high-density regions than in low-density regions, but this con- 
clusion is based mainly on the value of a single density bin, for 
overdensities from 1.8 to 10. The large size of this bin may mask 
the higher fraction of red galaxies in the upper part of the bin. 
Although, because of our much smaller galaxy sample, it is diffi- 
cult for us to measure the overdensity in exactly the same way as 
CNC07, we obtain a mean overdensity of >8 for the galaxies in 
CI 0332-2742. The evidence of a red sequence in CI 0332-2742 
seems therefore at variance with the conclusion by CNC07, but 
the degree of inconsistency is difficult to quantify. The disagree- 
ment may be explained if we are sampling a regime higher in 
density than studied by CNC07. 

5.3. Dependence of galaxy properties on galaxy density 

We have investigated the dependence of galaxy properties on 
galaxy density by comparing the ensemble properties of several 
samples (spike, field, high, and low density) and by relating the 
properties of individual z = 1.61 galaxies to the measured local 
density. The latter method did not provide significant evidence 
of correlations between the measured galaxy properties and the 
local galaxy density. There may be several reasons for this lack 
of significant correlations. First, the local galaxy density is de- 
termined from the secure spectroscopic members and the less 
secure photometric members on a scale of -300 (physical) kpc, 
which should be optimal for the detection of galaxy groups, but 
is rather small given the mean distance to the closest neighbour 
among the 42 confirmed members in the GMASS field, which is 
of the order of this scale. We note that, when comparing galaxy 
properties, only the spectroscopic members were considered, al- 
though the photometric members have a significant influence on 
the computation of the local galaxy density. For example, the 
strong density peak in the western part of the GMASS field con- 
sists almost entirely of photometric members, which certainly 
merit spectroscopic follow-up. It is also possible that the influ- 
ence of galaxy density is notable only on larger scales, such as 
those of the high density region, for which we do find a signif- 
icant correlation between density and all measured properties, 
most notably morphology. 

The environmental dependence of galaxy properties was 
studied at higher redshift for a significant overdensity of galax- 
ies at z = 2.300 by Steidel et al. (2005). They found galaxies in- 
side the overdensity to have mean stellar masses that were higher 
and inferred ages that were ~2 times older than identically UV- 
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selected galaxie s outside the struc ture. This is similar to our re- 
sults atz = 1 .61 . Peter et alJ (l2007l) augmented this analysis with 
one of z = 2.300 galaxy morphologies but did not find evidence 
of a diff'erence between the proto-cluster sample and a control 
(field) sample. Very few of the z = 2.300 galaxies appear to be 
regular ellipticals or spirals. The morphological diff'erence be- 
tween overdensity and field, which we found at z = 1.61, is due 
to the elliptical galaxies, which are also red and, using the UV 
selection applied at z = 2.3 00, could therefore not have been 

found by [SteideletalJ ('2005'). 

The lower SFR that Steidd 'etaD (l2005b and ourselves found 
for galaxies in the overdensities compared t o those in the field 
appears to contrast with the conclusion of "Elbaz etaLl (l2007l 
EDB07) that the average SFR derived from FIR and UV pho- 
tometry increases with galaxy density in the redshift range 0.8 < 
z< 1.2. However, the range of projected galaxy densities probed 
by EDB07 is between 0.4 and ~3 Mpc"^ in recessional velocity 
bins of V =~ 1500 km s"^. At z = 1.6, this corresponds almost 
exactly to the redshift bin 1.600 < z < 1.623, encompassing 
all 42 spike galaxies. Using this bin and the rather large boxes 
of 1.5x1.5 Mpc^ used by EDB07, the projected galaxy density 
in the GMASS field is 5 Mpc"^ on average and 10 Mpc"^ for 
the highest density regions. Obviously, using smaller boxes for 
counting the galaxy density, as also noted by EDB07, would re- 
sult in even higher projected galaxy densities (up to ~40 Mpc"^). 
This shows that we are probing higher densities than studied by 
EDB07 and their conclusions do not apply to CI 0332-2742. 
Although this is probably not a statistically significant result. 
Figs. 8 and 9 of EDB07 appear to indicate that the SFR does 
indeed decrease at projected galaxy densities of >3 Mpc"^. 

6. Summary and conclusions 

In this section, we summarize the results presented in this paper 
and add a concluding remark on the galaxy structure found. 

- We carried out a spectroscopic survey, called GMASS, 
in a part of the CDFS, targeting galaxies detected with 
IRAC/Spitzer at 4.5 yum and pre-selected on the basis of a 
photometric redshift z > 1.4. We obtain spectroscopic red- 
shifts z > 1 .4 for 135 out of the -200 observed sources (Kurk 
et al., in prep.). 

- The redshift distribution of these sources shows a prominent 
peak at z = 1.6 consisting of 42 galaxies. We verified that the 
redshift distribution is not significantly aff'ected by observa- 
tional biases and find that this spike represents an overden- 
sity in redshift space by a factor of at least eight. We call this 
structure CI 0332-2742 and the galaxies at z = 1.6 member 
galaxies (of this structure). 

- The angular distribution of the member galaxies is not ho- 
mogeneous: within the 7'x7' GMASS field, there is a region 
of about 3 Mpc^ where the number density of members is 
five times higher than in the rest of the field. We call this re- 
gion the high density region. The brightest member is located 
760 kpc away from the centre of this region, within 65 kpc of 
two other bright members. These galaxies are not only the 
brightest but also the reddest members. 

- From the 42 redshifts of the members, we determine a ve- 
locity dispersion of 440+^^ kms"^ for CI 0332-2742 and of 
500^ kms"^ for the 21 members in the high density re- 
gion. The redshift distribution of the 42 member galaxies is 
best fit by the sum of two Gaussian functions with a diff'er- 
ence of Sz = 0.008 (or 920 kms"^ in the reference frame of 
CI 0332-2742) in their central redshifts. Although we do not 



find evidence of a spatial separation between the two groups 
of galaxies centered at z = 1.602 and z = 1.610, we cannot 
exclude the existence of substructure along the line of sight. 

- Using various methods to derive the mass which will be or 
is contained within CI 0332-2742, we find a mass within the 
range 6 - 60 x 10^^ M©. The irregularity in angular and red- 
shift space, including at least two localized higher density 
peaks, suggests that the structure is not yet virialized but will 
evolve to become a galaxy cluster at a later cosmic time. 

- We used the 1 Ms Chandra observations to search of hot gas 
associated with the high density region but did not find ev- 
idence for an X-ray luminosity > 3.5 x 10^^ ergs"^. This 
X-ray luminosity, however, is consistent with the expected 
X-ray emission for a cluster or group of galaxies with a ve- 
locity dispersion of 500 km s"^ 

- In a colour-magnitude diagram, a bimodal distribution of 
member galaxies is seen with blue galaxies around z - / ~ 
0.6 and red galaxies, including all early-type members, at 
z - J > 1.3. The latter form a red sequence consistent with 
model galaxies composed of stellar populations formed at 
a redshift z = 3. The measured scatter around the red se- 
quence of 0.135 magnitudes is consistent with the expected 
scatter at this redshift. The bright tip of the red sequence is 
formed by the triplet of bright elliptical galaxies mentioned 
above. We speculate that the clear bimodality seen in colour 
for CI 0332-2742 members, which is stronger than in field 
samples at this redshift, is due to the bimodality develop- 
ing earlier in volumes of high galaxy density, consistent with 
galaxy formation models where massive galaxies form first 
in the strongest density peaks present in the early universe. 

- We compare various galaxy properties of the sample of 
member galaxies within the high density region with the 
sample of member galaxies outside this region and with a 
similarly sized field sample of galaxies at spectroscopic red- 
shifts determined by GMASS in the range 1.4 < z < 1.8, 
but excluding the member galaxies at z = 1.6. These galaxy 
properties include the galaxy type based on spectral fea- 
tures (either early or late), the galaxy morphology (either 
elliptical, spiral or irregular), broad band B - I colour, and 
four properties derived from broad band SED fitting: mass, 
star formation rate, specific star formation rate and age. The 
galaxies within the high density peak are on average older, 
more massive, and redder than those outside the redshift 
peak. The SFR and sSFR of the members in the high den- 
sity region are on average three and five times lower than 
in the field. These results are consistent with current views 
on galaxy formation and evolution, where the most mas- 
sive galaxies form preferentially in the highest density peaks 
of the cosmic galaxy distribution, which collapse first. The 
massive galaxies that we observe in the overdensity also 
seem to evolve more rapidly than their lower mass counter- 
parts: they have lower specific star formation rates, already 
at z = 1 .6, which implies that they have formed most of their 
stars by this epoch. 

- We construct composite spectra of blue and red, early- and 
late-type members and of members in the high and low den- 
sity regions. The composite spectra show features consistent 
with the results obtained above: the high density composite 
is clearly redder than the low density composite. In addition, 
the [O II] line equivalent width, a measure of star formation, 
is clearly lower in the high density composite than in the low 
density composite. 
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We note that this is the first and only structure of this kind 
known: its redshift of z = 1.6 is higher than that of any known 
galaxy cluster and the structure contains spectroscopically- 
confirmed red, early-type galaxies, whereas at higher redshifts, 
the known and possibly only members of galaxy overdensities 
are blue, star-forming galaxies are known (such as LBGs and 
LAEs). CI 0332-2742may therefore represent a link between the 
z < 1.4 galaxy clusters and z > 2 galaxy overdensities. The ex- 
tent of this structure is not limited by the size (-I'xl') of the 
field covered by spectroscopy carried out within the GMASS 
project and neither does it see m to be limited by the size of the 
CDFS (Castella no et al.l l2007). Its extent of more than 10 Mpc 
therefore suggests it is a sheet in the web-like distribution of 
galaxies. The region with the highest surface density of galax- 
ies within this sheet, called CI 0332-2742 in this paper, already 
contains seven massive, passively evolving galaxies and already 
has a velocity distribution typical of a large group of galaxies in 
the local Universe. Since the number of passive galaxies, its ve- 
locity dispersion, and its mass can only increase during the -10 
Gyr between z = 1.6 and z = 0, reaching properties typical of a 
present-day cluster of galaxies, we conclude that, by observing 
CI 0332-2742, we are witnessing the assembly of a cluster of 
galaxies. 

Appendix A: Spectra and images of the member 
galaxies 

Figures lA.ll IA.2I IA.3I and IA.4I show the spectra and images of 
the member galaxies. 
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Fig. A.l. Spectra and postage stamp images of the 42 member galaxies. Wavelength in Angstrom on the horizontal axis and flux 
in 10"^^ ergs"^ A"^ cm"^. Uncertainties caused by background noise are indicated by the underlying filled grey spectra. As the 
flux calibration of the GOODS spectra is not very accurate ( Vanzella et al. 2006) but the absolute flux level of these spectra seem 
systematically higher than the GMASS spectra (for the same magnitu des), they have been divided by a factor two. The origin of 
the spectra is plotted: GMASS (Kurk et al., in prep.) or GOODS (V anzeUa et al.l(2006l) . the GMASS identification numbers, and the 
spectroscopic classes (see text). Robust (tentative) spectral features are shown by solid (dashed) lines with numbers (see below). 
The postage stamps are constructed from the HST/ACS observations in the B, V, and / bands, convolved with a Gaussian kernel. 
Indicated are the morphological class (on the left) and the z magnitude (on the right). The spectra are sorted on z - J colour, 
starting with the reddest galaxies. The numbers for the spectral features refer to: (1) Si IU1527, (2) C IV/11550, (3) Fe 11/11608, (4) 
Al 11.11671, (5) A1IIU1855, (6) CIII]^1909, (7) FeIIU1926, (8) CII]^2326, (9) Feim344, (10) Feim375, (11) FeIU2383, (12) 
Feim587, (13) Feim600, (14) Mgim796, (15) MgIR2804, (16) Mgm853, (17) [OII]^3727. 
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